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ABSTRACT 

This work describes a method of determining the ratio of heat anda 
momentum diffusivities from micrometeorological data obtained by the 
Great Plains Turbulence Field Program bAR 

The method involves obtaining first the value of the Deacon sind 
Meer le parameter, usede&s the indicator of stability, and with this 
value, the ratio of two forms of the Richardson number, which is theo- 
merteaiily equivalent to the ratio of heat and momentum diffusivities. 

This investigation was carried out at the U. S. Naval Postgrauuate 
Semool, Monterey, California, durimg the period February 1959 - Hay 
1959, in partial fulfillment of the requirements for the degree of 
Meter of Science in Meteorology. 

The writer wishes to express his appreciation for the assistance 
miemezordance given him by Professor Frank L. Martin in this investiga- 
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ie lene POGe LON 

A cneory, 1s ieeetcedteancer Lettau [6), and Martin (aur 
which permits computation of special forms of the Richardson number 
and flux Richardson number. These equations are dependent upon 
the Deacon wind profile parameter C > and the roughness para- 
meter Z2,. Equations for the computation of these parameters 
are derived also. 

A recent theory is included in these determinations, and in 
a sense, this paper constitutes an attempt at testing this 
theory and comparing its results to that of the many previous 
investigators in Che fielavel atmospheric turbulence. 

The data utilized in this study is evaluated and is considered 
to overcome many of the considerable limitations inherent in data 
evalileple {or earlier research, 

Correlation coefficients of the resulting ratio of heat and 
momentum diffusivities with the two forms of the Richardson 
number and the Deacon wind profile parameter were computec,and 
the dependence of the ratio upon the other variables is discussed 
and compared in the light of knowledge gained by earlier investi- 


gators. 





2. Theory and derivations 


The work of von Karman, Prandtl and Nikuradse has shown 
that the fundamental relationship for the wind shear within a 


neutral surface layer is given by (2], 


34: Huo. (1) 


= 
= 
= R 2 


where u = the mean speed at a height z above the surface 
uy, = Vex called the "friction velocity" in the neutral 


layer 


se 
iH 


the shearing stress at height z 


fluid density 


p 


k von Karman's constant 


Il 


On integration, Eq. (1) gives the logarithmic law for the 
neutral surface layer (defined as the layer in which the shearing 


stress is effectively constant), 
Ma = Age In 2 2) 
e 


which applies to fully rough flow, and 2, the roughness parameter, 
is the measure of surface roughness. 
In studying the wind profiles versus the logarithm of 
height as ordinate under non-neutral conditions, Deacon (2) 
noted that the wind profile differed from the logarithmic profile 
in that it showed convex curvature under unstable conditions, 
and reverse curvature under stable conditions. He has suggested 
Cnae tne, Telationcmip, g : 
$52 fos * (3) 
Rai" 
should exist under non-neutral conditions, where @ is the wind 


profile curvature parameter. 





Applying the condition that u = 0 at z = fos Which hoids 


for the logarithmic profile, and integrating, this expression 





becomes, | 
Mx (( 2) Ler) (iL | 

AR(i- BY Zo y 
Assuming that kx | remains constant in the layer z, to 


R= 9) 
Zo, one may write (2), 


U2 (2/2) °- 1) f= 
at > eT (5) 
Determination of 2, using an atmosphere of neutral stability 
allows one to calculate @ knowing the winds at the heights 1 
emda 2. ha. eee as known, may be solved for 6 
( 
by a method of successive approximations known as the Newton- 


Raphson method. In this method, successive approximations to a 


root of any equation f(y)=0 may be written in iterative form as, 


Hie Sp Se ( ) 
int — Fu TF Q 
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Mearranetliguags (>) Tnto cue form iy j=05, were. y-i. = @ , one 


obtains, 
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Placing the above expression into the iterative form, and making 
e of 2 ‘2 ol 
msc OL (the tdentirty, Deg. on ios bee (S ) Legho 


sives, - eer ‘i _ (J+ 17 
‘ Ze ~ 


This method of obtaining es was used in this paper and 
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requires only several iterations to obtain the accuracy specified 
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Since the solution for @ from Eq. (8) is uependent 
upon the roughness parameter, a method of determining its 
Value under the applicable conditions must be deyiced, As 
inuicatea, 2, values may be cetermined under neutral condi- 


byOnG a Fount Se fica Stnse.. Since ie 


‘ MAdCPCHAeNL “OF 


pts 


lapse rate [3] , use is made of the following cinmensionless 


ratio, due to Davidson and Barad [1] 


> 


 - Se Boe, 2 es (4) 
[yr Uz alec jp ~ SEU 





lato Wiel Ene expression. for Cue logaritnmic: jroli teem, yc 
substituted. This gives 

) f f 0 

aul —— an 2,tdn2: - AAjtes ~ 7 Az, {O 


ZU Indy +dn2; thin 2, + tn 2, — Pn 2o 
Gpon Lurther reduction of Eq: (10) an expression :62 7258s 
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where Yn2; 2 dye dn2gtdn2 tne, 5 Ly Mrz In 2,tdin2s nace, 


de, (urbe lenec, Eheoary 

Ie CGELVatiOn Of ENG Penainde, or Ehevenegry sin skis 
paper is dependent upon Lettau's turbulence model [6] , 
Walch Was adapted to the data in a manner similar to that 
of Martin [8]. 

In avatogy to moicewtar cifiusion, tEhe theory ofstwbu- 
lent mixing defines the exchange cocfficient as, 

A=~ @hiw (1 2) 

which is written by Lettau [6] as, 








ate 
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where w is c#tts) the mixing velocity and 1 the root-mean-square 
mixing length introduced by Prandtl {3}. Since turbulence 


theom also states that, 


_— 


Tz - guts AC%z) (14) 





one may write that, 
wT : ' 
f= et we (15) . 
we 
meere Ww 1s called the virtual friction velocity. 
Combining the last three equations we obtain, 
6a yu io 
&= L( 02) (16) 
Lettau (6] assumes under adiabatic conditions, that 


Ux Be acer (7) 


am Cg 


Under non-adiabatic conditions, Lettau [6) makes the following 


assumptions, si oe 


a (W8) 


co ae 


we & x” Z ae 

| ete ann We _ N Cc C CG ( q . 
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where Ng= Kyi Rays the ratio of the coefficient of eddy diifusi- 


mrty for heat, to that of eddy ditiusivipy Lor momentum. 


ats 
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Substituting the expression for wv from Eq. (18) into 


Eq. (19), gives upon rearrangement, 


) - te _ Nghe 26 (20) 
X 


in 





Peeshould be noted at, this point, thatermethe last term of 


Eq. (20), the term Ng represents the efficiency in the conversion 
Of “ermal cnergy into kinetic energy. The last term of Eq. (20), 
which will be denoted x, may be written 


x =NoX, X= gece 2 (2.1) 


Cm 


NY. 


where X is a particular form of the Richardson number introduced 


by Lettau [oJ]. Actually, 


whereas 





Lut X will henceforth be called the "Richardson Number", 
x will be called the flux Richarcson number [6]. 


Combining Eqs. (18), and (19), Lettau (6) obtains, 


& - fe (? 3) 
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Gombining tose (3), Gla 46 )emand (23) enab lesceone gic 


Gegive the followin results, 
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x t-2P 
de Wo (2 (2 5 ) 
This derivation is shown by Martin (8). 
(bs “Use sot finite-diirercnce Cerivatives 


In computations with Eq. (25), Lt is necessary to wee 


fanite-difttercnee derivatives. The instruments were arranced 


at constant logarithmic increments so that the following 





transformations due to Lettau [7J arc of value, 
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n> = 22/3, (Z 6) 
2) 2> crea, te 


where 4“, 9 represents the geometric mean of any layer yy LO 22, 


and nié, for the data levels utilized in this work. Therefore, 
2 ( 


ae Pn tCrms Of finiveqrl ferences. 9445 becomes, 
ee 5 Up-UW. Uz-U- -% (2. &) 
O2 2)n 2s z Lin z lind 


Letting 4=21 5 makes Eq. (28) representative at a mean height 


in the layer z, to 2), and defy, becomes, 


Un- U ~ We n\n" (2 q) 
Z2Imh ~ fea (lex) ae 


Rearranging leads to 


pet > Fede (an)! G0) 


re) > (Ua)Rk Va, 


and solving for x, gives, 


=| "al . 2 ) 
Vv nd (= 2 - =. i 
(iu ) 
“< 
Hence, it becomes possible to solve for x, knowing Wi, @, and Zo. 
(ce) An expression for wa 
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Sicictly speakin® Ve Var mesa sta mere y anc 71 bu ceo 
a ref » 
Stroplic wind. \fhe Variation Wren stability is-notlye. accu- 


rately known put the variation from one adiabatic time to 





; , : “ ; 
another may be assumed to be linear. An expression for Le 
terms of the neutral-wind profile may be obtained by lovarithmic 


differentiation of Eq. (2), with Ura replaced by wa, giving, 


Pa a oe Cn ee, Va 
Oo kn a oa ‘© J) 


whici son expansion, becomes 


rr 


a ee | Uay — Wai) + (Ug, - - Usd) _ 
, Gna, “tna, JF Ome, =e, )n ame _ 3) 


where (Up. is notation for 
. Dae: oe ‘he 
Ni Ug, + Ua, ~ Ua, = Aa 
ane yaa , 
lopli lo on = 41n2 with instrumentation as set up in the 
Greate: ialiis project . 
Substitution of Eq. (33) into nq « (31 gives 


ees = {| Aes so} (24) me i} (3 "2 


Variation of Au, linearly to coincide with the bi-hourly observa- 





tions between’ the neuenal times at sunrise and sunset will give 
representative values of the™Elux Richardson number 3% during 
the day. 
(d) <A convenient formula for X 

The particular form of the Richardson number used in the 
Lettau theory was shown in Eq. (21). This form, hence forth 


denoted X, may be written as [7], 
é 
in Lo. 


Similarly to Eqs. (26), we have the transformation 


ge We . te 
20> 25>: Se (3 6) 





end tiuerefore in terms of finite-differenees, the followiitz rela- 





ai" _ Cee © : on _ fat - Ce w{ ie T) 


2,.,4mn ? : a, oe, n* 


SUeoCilLucinegehds. (3/7) into Eq. (25) @ives, 


? 2 2,.(0,-@,)ln n> 38 
x= 5, (Ua, -Ute,)* G ) 


memend a5 = K:,, may be eee over all successive pairs of 
levels to be A4A4 of Eq. (33). In Eq. (38) Om may be approxi- 
mated by Oy = T2. Making the substitution lk? - Ua = AuU,/% gives, 


— (bg WnAlO,-C,) 2,2 =. Ch 
o> 8a G4) 


VW. 
Siace hid no longer appears in either x or K, but has been 


~~ 


qe 


replaced by AU,, the latter has been varied linearly between 
times of neutral stability (near sunset or sunrise). 
(ec) Ng as a function of height and other variables 
In Eq. (21), it was shown that, 
Niw= =e 
Cu” x 

and we now have convenient working equations for x and X, 
Eqs. (34), and (39), respectively. Computations of Ng as a 
funetion of elevation can now be made provided the roughness 


parameter Zg and the wind profile parameter @ have first been 


obtained. 





Pees eure Of CN Gata 


In addition to testing the micrometeorological theory set 


=F 
lar 
Nd 
th 


cth in this paper, this tavestigation was also initiated by the 
availaLility of comprehensive round-the-clock data of the nighesi 
caliner to date. The Great Plains Turbulence Program [£7] was 
a research project of the U. S. Air Force Geophysics Research 
Directorate conducted in open prairie country near O'Neill,. 
Nebraska between, 1 August, and iQ September 1953. 
POUERtLCLEN Organizations participated in thewpeajcct .aanic 
their contributions included: -instrumentation, data collection 
and evaluation, development of turbuience theory, instrument 
description and evaluation, site preparation, and data tabulation. 
The Great Plains Turbulence Program was in progress ior about 
B1lx weeks Gering which time seven observation periods vere com- 
begsteG. The criteria for an acceptable observation periou were 


trat the ca,time periods be clear to partly cloudy, with clear 


oe 


skies at night, so that large diurmal variations in the heat budget 
terms, and particularly, in the Richardson number could 
Also, it was preferred that the megan geostrophic wince speed vary 
from one observation period to the next. Instrument arrangement 
on the £ield site required that the general observation periods 
be charaetcrized by relativels constant wine €ireéction betveen 
SE and SW [7] . 

Only the data of the sixth and sevemeh general opeeraimo: 


parigds, for the days, 3h August through | September, and 1-3 


Sectember 1953 CY) , respectively, will be described anc cvaluatec 


10 





in thie paper. 

paesaciuoeers Institute of Technology tad Joms Heokins 
Uitvemertty (Weredmmer tClérrea to as Mit, and JHU, G@escetively ), 
wind and temperature profile observations at standard Levels were 
uUbplized, The Mildata collectedvarethe Zz. 4, ServandwiG-meter 
levels consisted of L5eminute mean wind speed (cem/see}, and 157 
minute mean air temperatures (°C). JNU data was collected at 
, »6, 1.6, and 3.2-meter levels and consisted of Wourly-meaa 
wind speed (cm/sec), and average temperatures From 20 readings 
@uriny a S7~minute intcerWwal., Bi-hourly obs@rvations were taken 
from the. time the general observation period began with ali tine 
averages centered at observation Cime. 

The standard cup anemometers used by MIT were of tne conven- 
tional type ana vere wind-tunnel calibrated prior to installation. 
bifferences in calibration were minimized by field matchine tests 
cundWetea at the test site. Wind measurements were averazca 
over loeninute intervals and recorded to the nearest em/sec. 

JHU obtained a number of three cup anemometers of 1341-1942 model 
from the Soil Conservation Service, and by modification and careful 
matching produced a set of mechanical wind measuring instruments 
with compatability of the order of 0.3%. Wind measurcmente were 

of hourly mean wind speed reeorded to the nearest cm/sec. 

MIT temperaturcs were ncasurcd with identical, artificially 
ventilated thermocouples, similar instrumentation having previously 
been usce et MIT's Round Hill Field Station. The thermojunctions 
reve shielded from sular and ground radiation and ventilateu at 


ayer’ ok > mp Sec. Calibrations CONGLCLCa 11) Wace, patiis 


- 
tl 





matched against precision mercury thermometers. Fiftceneminute 
mean temperatures at the standard levels used in this paper were 
ree@mued to the nearest 0.01°C. 

The description of JHU temperature installation is similar 
to that of MIT's, the main differences being in the method of 
shielding awd the lack of artificial ventilation. Calibration 
was similar to MIT procedure previously described. Using Beckmann 
differential mercury thermometers between the two water baths, 
aieabsolute accuracy of 0.05°C was indicated. An average of 20 
temperature readings at each level during a Seminute period was 
reecordcd tosthesmecarest O,OI°C, A detailed description, of thu 


iw@etrumentation, ectc., is given in Lettau: (73. 
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G4. Jomputations 

The solution of the working equations for @ , 29, X, and X; 
equations (8), (11), (34), and (39) respectively, was accom- 
Peisied by desk caeulator, rather than*by slide rule or eclectro- 
nic computer. It is felt that use of the desk calculator was 
justified considering the relatively smal quant liyeor data. 
Also, greater accuracy than that attainable by slide rule, was 
aosuved, to ttWe limits of the@@ata signi ftteance. 
(a) The roughness parameter, zo 

Determination of 29 by means of a fitted straight line on a 
in 2 versus Wind speed graph can afford only @ rough approxima- 


@ifen, at best. Aecordingly, the worl@ime equation, 


= Sues Snr —_vr 
iL 4 Eu 
applicable at the adiabatic times of sunrise and sunset, was used. 
Fortunately, the 0635 and 1835 observation times were near enough 
to local sunrise and sunset at this time of the year (31 Aug. - 
8 Sept. 1953) to satisfy the adiabatic requirement. 


For the MIT data levels of 2, 4, 8, and 16 m, Eq. (11) with 


24 = 24, 23 = 22 AO} 22 | and z}] = 2 meters may be reduced to 


2. = bee \(2.5 = Ee)leg 2 (40) 


Tie roughness parameter for EhéevjHU data levels of .4, 9.6. 176, 


the form 


dndes.2 Mm, 1S.obtained by reducing shq, (1). ta. 
an i; ach uh 2 41) 
ae ney Leg, (0.1780 1S ar zu ) hegre A 


where avain, 2Us Uyer Us + 4+, and, AU= UW, Uys As - OY 
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The above formulas are solved very simply by desk calculator 
and logarithmic tables. The value of Z, obtained by this method 
is made to hold for all computations within a 12 hour period 
centered on the time for which it is determined. 

(b>) The@eaeerativye procedure for yi = 1 = @ 

Eq. (&) required only that the roughness parameter be known, 
in addition to the observed data. No suitable objective method 
fom obtaining an initial estimate of y; could bem@evised. Accord- 
ingly, after some experience in the use of iterative steps, it 
was found that subjectively-determined starting values for the 
other bi-hourly times, based on the premise that @ increascs 
with diurnal heating and approaches one nearer the ground, provided 
for more rapid convergence of the computation. It was observed 
that the solution converged quite nicely for all but near-neutral 
conagitions, where considerable difficulty was encountcred. 
mecording ly, iriteria for accuracy of the solution weme devised 
which simplified the computations a great deal. That is, if 
lg sh? sO , then the solution could be halted when two places 
of decimals had been fixed. However, for Iguhs .t0 , and if 
the solution was oscillating over the desired root, it was 
required that three places of decimals be unchanged after further 
iterants are obtained. If the solution was not oscillating, 
but approaching the root from above or below, then four places 
of decimals were determined without change in a successive 
iterant. Also, it was decided to disregard the solution if 


ly31¢ O05, for this case was considered identical to @¢@ =l. 
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A rary 3 to A eS oe 
steluteak SCBE so Gee Ba 


In cadition, certain uata-series of hicar 


Were exeluccu: in particular, those case; tu which 46,-@,| 


in ny layer was less than 0.05°C. 


(cc) Conputaion of x 


As Was pointed out in the discussion fullowing Eq. (31) 


the solution ror x reqtires the knoWhedsem? 7), e@ 


rd 
representation of wg under non-neutral conditions, This repre- 


2entation vas accomplished in Eq. (34) by allowing the function 


Aus Uy 7 Us. ) + (Un, - Uc, ) to vary linearly between tne 


Peer 2nd losgoeadiabatic conditions. The vaive of % cl meincd 


is generally netative for unstable cases and pusitive for ptable 


The roughucss parameter has already been ceternmined for the 


wvilution of @ in any layer, and is applied in the solution for 


MeivotNe Sane manner as: betore.. The Vvalucwotee omar eet tae 


UrCr a, Li Ly. (J+) makes Enc value.of = Obtsinedvan licen! a6 


foc Meontlric mean of the laye. concermed, wiiercas tlic @ value 


4 


deol CO DC FeprCcentative Ol Ticwewbr ree ere) veoee a. 


WAS 
oor suLatian ol x 


Tiesoluvion. of EG. (39) for the Richardson number cior any 


Layer Ls quite straightforward. As observation shows, it 


contain. only constants and values already ootainec, or 
fie cata. Tne value for g was taken to be JGO cm/sec”. Tc 

Ns a i : ' re) eS =e 
PON ao Vell COG M00 a2) Oi Bie expr cus orm > een Mele 


ey eS 
Cty ic ) ime Coos, ) Fi Was nsec where ; Bh, eae 


AS was mentioned earlier, a.criterion tor determining a non- 





neutral layer from the temperature sounding was set up, namely, 
{o,- ©, | aon Cin Che contrary case, the data was not used 
for computation purposes. The value of X obtained has the same 
sign as the temperature lapse, negative for unstable conditions, 
MmienpoOssinvembor stable conditions. Since 2] > appears in the 
equation for X, its solution is also applicable at the geometric 


mean of the layer. The working equation for X now becomes, 


sp - 1b gt all(T.-T) +2) Me. U2) 
Te (AU) 


{ 


Fhe Sotution for Ng is Simply Ehe Vvalvicrwok x iromiikda., 1¢579 
divided by the value of X from Eq. (39). The value for Ng 
thus obtained is always positive, and applies at the veometric 


Wea Of “the Vaver, 2] to 29. 
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Seer SCUSS LON On Lesuits 


This section contains a summary of 


results for the computa- 
tions of g SH and Ng» and compares these results to those 
obtained by previous investigators. The numerical values of the 
avove parameters are shown in Tables 1, 2, and 3 for the indicated 
times, and layers. 

Deacon (2) quotes values of the wind profile parameter, which 
he obtained from Eq. (5) over short grass, using wind speeds for 
heights between the 0.5 and 4emlevels, as ranging from 0.72 to 
Pets. Ymspection of Tables’ 1,2, and 3 shows a range of @ of 0.865 
to 1.14, a somewhat lesser range than that obtained by Deacon [2] 
trom the series of data he utilized at Porton, England. 

Priestley [9], shows a plot of log Ng versus x, for mainly 
unstable periods, in which values of Rf range from 0.1 to -0.5. 


His scatter diayram is shown below as Fig. 1 of this paper. 


Based on this data, Priestley obtained an average correla- 
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Fix. 1. Dependence of Ky/KM on Rf at 1.5 meters 
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Computations from the sixth general observation period 
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Table 2. Computations from the seventh general observation period 
(MIT data only). 


Time Layer X X N 
ee € 9 
Wess 2-4 Osa = 5)2)50, - 0924 iS 
1435 8-16 1.1041 - , S013 = Oo RG 9.669 
a 4-8 1.0957 -.9463 —wOoIG Yo TOL9 
bs 2-4 ie 0050 -,9160 - 0924 9.909 


Table 3. Computations from the sixth general observation period 
(JHU data only). 


Gime Layer GC x X Ng 
(m) 

1035 Oso 2 1.9394 -.2192 =. US. Blige pels) 

1235 3 L.0657 =. Uos 0005 6.262 

1035 ~O-1.6 1.0299 =.124/ = OSU) 2.200 

1435 1 oe errs 1.0490 a a 12) - .0683 4.365 


Table 4. Computations of the roughness parameter 











Date Time Zgo (m) (MIT data only) 
31 August 0635 0.017960 
. 1835 0.012537 
1 September 0635 0.015483 
7 September 0635 0.174535 
: 1835 0.013689 
Date Time 4%, (m) (JHU data only) 
31 August 0635 0.008249 
. 1835 0.007335 


a a a a a eaten 
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A similar comparison made from the results of this investi- 
gation, shown in Figs. 2 - 6, indicate that the range of Nog is 
several times as great as that shown by Priestley [9J, or as 
calculated by Lettau (7J from the Great Plains data. Likewise, 
values of x computed in this paper ranged up to af 1.0. Exam- 
ination of Priestley's scatter diagram, Fig. 1, shows that his 
range of x <0 is much smaller (closer to neutral conditions ) 
than the points for which x<0 in Figs. 2 - 6. There is also an 
indication in Fig. 1 that his curve is making a bend in the 
vicinity of x=0, similar to that on the right side of Figs. 2, 
Be and 4. 

It was noticed by this investigator that the computation 
of x was quite sensitive to erroneous values of the roughness 
parameter. In the case of very rough surfaces, the vertical 


wind profile may be represented in the form 


~ Vs ln 2-4 
— R a5 


iA 
where d is called the datum~level displacement (3). Lettau C7] 
makes use of the zero-point displacement D=d + Zo in computing 
values of Z at the JHU site, and obtains an average value of 
d=9.9 cm. This means that the logarithmic profile produces a 
zero value of wind 9.9 cm above the zero point of nominal height, 
Zo. In the case of very dense vegetation, d+ 2, 1s equal to 
the height of the vegetation. Lettau C7] reports the value of 


d he obtains to be on the order of the grass height, 10 cm. 


Lake [5] , and Deacon [2] also make use of D in computing values 
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- Ng versus x for MIT aata in the 8 - l6-m layer. 
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Ng versus x for MIT data in the 2 - 4-m layer. 


Figure 5. 
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Figure 6. 





of Zo. Deacon gives values of zp for mown grass=surfaces of 
grass length 4.5 cm, ranging from zg = 2.4 cm to ¢g © 1.7 «em, 
depending upon the strength of the wind at the 2em level. The 
results of this investigation gives values of z, ranging from 
0.73cm to 1.80cm, as shown in Table 4. 

Lettau [7j calculated values of X uSing an equation similar 
to Eq. (38) at the 2,4, and 8-m levels. Values of X obtained 
in this work are consistent with his results and are valid at the 
geometric mean of the layers between 2, 4, 8, and 16 m. 


Results of correlations of N@ with x, X, and @ , in the 


entire layer 2 m to 16 m are shown in Table 5. 


Table 5. Correlations with Ng with x, X, 


and G : 











Sample Correlation Coefficients Significance 
size Ng vs. x Ne vs X Ne vs G Level 
Stable 10 965% »/07* o911* 939 
Unstable 21 -.805* aoe =, 244 = 440 





Note: Significant correlations for the 0.05 levei of belief are 
indicated by the asterisk (*). 


The 5% significance level was derived from the test men- 
tioned on pp. 124 of {4] . All of the correlations obtained 
under stable conditions were positive and significant, indicating 
that Ne varies with stability in the same way with x, %, and o 
However, the correlation was highest with x, which is to be 
expected since x explicitly contains Ng 1n its formula. In un- 


stable layers, all of the correlation coefficients were negative 


Q 


o~! 


although only one of these was significant at the 5% level. 
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This particular correlation verifies,in a general way,that of 


, 
Priestley [9] ,also taken under unstable conditions. 
Fig. 7 shows a plot of Ng versus time for each layer in 


which N@ was determined. For the period 31 August - 1 September, 


the average values of Ng are listed below. 


Layer (m) Ng ave 
oe eo ibe Leo od 6 
4 - 8 4.9735 
Mn Bei D530 


Note that Priestley's (9) mean value of Ng at 1.5 m would be near- 


ly one. However, he has values of Ng as high as 4.0. 
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in the layers 8-16, 4-8, and 2-4 m, versus time, for the rerioc 
1 September 1553. 


PLOt2oFf Ng 
1 August - 
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6. Conclusions 

The solutions obtained in this investigation are subject 
to rather stringent assumptions which undoubtedly account for 
some variance of the results from those of previous investigators. 
It has been shown in Section 5 that values of Zoe 8 and X, Lie 
in the range of previous investigations, but Ne is somewhat 
larger than comparable values previously owtained. 

The theory presented in this paper is developed to allow 
computations of the various parameters for non-adiabatic condi- 
tions, since the ‘accuracy criteria of Section 4 ruled out ail 
computations for near-neutral conditions. By Eq. (34), errors 
in MAU, which tend to increase x tend also at the same time to 
decrease X [from Eq. (39)] , and therefore increase Ng: 

However, this investigation did not attempt to verify the 
above observation, but to put forth a method of computing the 
pertinent turbulence parameters. The technique of representing 

o/p by linear variation of the function AW, is a previously 
untested one and may be the contributing Pition to the varia- 
tions of Ng thus obtained. 

Since these computations ignored the use of @ zero.point 
displacement in computing values of the roughness parameter, 
no conclusions as to the contribution of Z, to values of Ng can 
be made. However, reasonable solutions for 2 were cbtained., 

The plot of Ng versus x in Figs. 2 - 6 shows that Ng 
increases, both with degree of ee and with degree of 


instability, and is near one for near neutral conditions. On 
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the basis of this paper it seems likely Fig. i would aiso show 
that Ng increases with stability if more points to the right of 
the origin had been obtained. 

The correlation of No and x showed a higher degree cf sig- 
nificance than that of Ng versus X. This was to be expected as 
Ng is explicitly contained in x. (see Eq. Ci Ae PriespLey 

{9] also reports that spuriously low correlations of Ng and X 
can be expected as errors in measured gradients affect Ng and 
X in an opposite sense. 

Consistent with the approach of this investigation, Section 
4 criteria ruled out solutions for many near+neutral layers. 
This investigation revealed that the seventh observation period 
(MIT data), and JHU data for the sixth observation period indicated 
neareneutral conditions much of the time; also solutions were 
unobtainable in the first five observation periods. Therefore, 
it can be concluded erOnmCNTS investigation that only the sixth 
observation period of the Great Plains Project contains data 
giving a significant range in ean ey, for the purposes of 
this research. 

An inspection of Fig. 7 shows that Ng is largest in the highest 
layer throughout most of the day. Secondly, Ng increases with 
increasing departure from adiabatic conditions in either sense. 
This accounts for the JHU data giving fewer results at its lower 


observation levels. 
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